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ground for nuclear physics, quantum chromodynamics (QCD), and nuclear superfluidity. Deter¬ 
mination of the fundamental interactions that govern matter under such extreme conditions is one 
of the major unsolved problems of modern physics, and - since it is impossible to replicate these 
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constraints come from measurements of NS bulk properties: each model for the microscopic be¬ 
haviour of matter predicts a specific density-pressure relation (its ‘Equation of state’, EOS). This 
generates a unique mass-radius relation which predicts a characteristic radius for a large range of 
masses and a maximum mass above which NS collapse to black holes. It also uniquely predicts 
other bulk quantities, like maximum spin frequency and moment of inertia. The SKA, in Phase 
1 and particularly in Phase 2 will, thanks to the exquisite timing precision enabled by its raw 
sensitivity, and surveys that dramatically increase the number of sources: 1) Provide many more 
precise NS mass measurements (high mass NS measurements are particularly important for ruling 
out EOS models); 2) Allow the measurement of the NS moment of inertia in highly relativistic 
binaries such as the Double Pulsar; 3) Greatly increase the number of fast-spinning NS, with the 
potential discovery of spin frequencies above those allowed by some EOS models; 4) Improve 
our knowledge of new classes of binary pulsars such as black widows and redbacks (which may 
be massive as a class) through sensitive broad-band radio observations; and 5) Improve our un¬ 
derstanding of dense matter superfluidity and the state of matter in the interior through the study 
of rotational glitches, provided that an ad-hoc campaign is developed. 
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Figure 1: Schematic structure of a NS. The outer layer is a solid ionic crust supported by electron degeneracy 
pressure. Neutrons begin to leak out of nuclei at densities ^ 4 x 10** g/cm^ (the neutron drip line, which 
separates inner and outer crust), where neutron degeneracy also starts to play a role. At densities ^2x 10*^ 
g/cm^, the crust-core boundary, nuclei dissolve completely. In the core, densities may reach up to ten times 
the nuclear saturation density = 2.8 x 10*^ glcrc? (the density in normal atomic nuclei). 


1. Introduction 

Neutron stars (NS) are the densest objects in the Universe. Composition evolves from ions 
embedded in a sea of degenerate electrons (and eventually neutrons) in the solid crust, through the 
crust-core transition where nuclei dissolve to form neutron-rich nucleonic matter, to the supranu¬ 
clear densities of the core, where theory suggests the presence of exotic non-nucleonic matter 
(Figure 1). The nature of matter in such extreme conditions is one of the great unsolved prob¬ 
lems in modern science, making NS unparalleled laboratories for nuclear physics and quantum 
chromodynamics (QCD). NS also host super-strong internal magnetic fields, and share the rich 
phenomenology of low-temperature physics. Although born in the hot furnace of a supernova col¬ 
lapse, NS cool rapidly to temperatures far below the relevant Fermi temperature. When matter 
cools there are two options: it can freeze to form a solid, or become superfluid (or superconducting 
if there are free charge carriers involved). NS are expected to exhibit both these phases. The outer 
kilometre or so of the star forms a solid crustal lattice of increasingly neutron-rich nuclei, while the 
star’s core hosts a neutron superfluid coupled to a proton superconductor and a free gas of electrons 
(as required to make the conglomerate charge neutral). The free neutrons permeating the crustal 
lattice in the inner crust are also expected to be superfluid. 

The most fundamental macroscopic diagnostic of dense matter interactions is the pressure- 
density-temperature relation of bulk matter, the equation of state (EOS). The EOS can be used to 
infer key aspects of the microphysics, such as the nature of the thr'ee-nucleon interaction or the 
presence of free quarks at high densities. Determining the EOS of supranuclear density matter is 
therefore of major importance to fundamental physics. However it is also critical to astrophysics, to 
understand not only NS but also NS/NS and NS/Black Hole (BH) mergers, prime sources of grav- 
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itational waves and the likely engines of short gamma-ray bursts (Bauswein et al. 2012; Lackey et 
al. 2012). The EOS affects merger dynamics, BH formation timescales, the precise gravitational 
wave and neutrino signals, any associated mass loss and r-process nucleosynthesis, and the atten¬ 
dant gamma-ray bursts and optical flashes (Metzger et al. 2010; Hotokezaka et al. 2011). The cold 
EOS probed with NS, and the way it joins up with the hot EOS that determines explosion condi¬ 
tions (Janka et al. 2007), are also vital to understanding the late stages of core collapse supemovae, 
including their gravitational wave signal. 

1.1 Unknowns in dense matter physics 

The properties of NS, like those of atomic nuclei, depend crucially on the interactions between 
protons and neutrons (nucleons), which are a starting point for ab-initio calculations of the nuclear 
many-body problem. While two-body interactions are well constrained by experiment, three- and 
more-body forces are a frontier in nuclear physics. At low energies, effective field theories based on 
QCD symmetries provide a systematic expansion of nuclear forces, which predict two- and more- 
nucleon interactions. There are also complementary efforts using lattice approaches to nuclear 
forces to provide few-body nucleon-nucleon and more generally baryon-baryon interactions, but 
this approach is still affected by large uncertainties. Some of the predictions of these models 
can be tested against current nuclear data; where data are not yet available, predictions are based 
on the consistency of the approach. The appearance of shell closure in neutron-rich isotopes, as 
well as the position of the neutron drip-line is found to be very sensitive to the three-body forces. 
Exotic neutron-rich nuclei, the focus of present and upcoming laboratory experiments, also provide 
interesting constraints on effective interactions for many-body systems. 

NS observations, by contrast, challenge many-nucleon interactions at extremes of density, neu¬ 
tron richness, and baryonic content. Astrophysical inputs are essential to make progress (Eigure 2), 
although extracting information from the data is harder than for lab experiments. CuiTent approxi¬ 
mate models all have major uncertainties, both at high density and for neutron-rich matter (Wiringa 
et al. 1995; Pieper et al. 2001; Stone & Reinhard 2007; Hebeler & Schwenk 2010; Steiner & Gan- 
dolfi 2012). Although experimental information on matter near the nuclear saturation density is 
plentiful, there are only a few experimental constraints at high densities and/or for neutron-rich 
matter (Tsang et al. 2009; Kortelainen et al. 2010; Roca-Maza et al. 2011; Piekarewicz et al. 2012). 
NS thus provide a unique environment to test models. 

At the very highest densities reached in NS cores, we also expect transitions to non-nucleonic 
states of matter (Glendenning 2000). Some of the most exciting possibilities involve strange 
quarks: unlike heavy ion collision experiments, which always produce very short-lived and hot, 
dense states, the stable gravitationally confined environment of a NS permits slow-acting weak 
interactions that can form matter with a high net strangeness. Possibilities include the formation 
of hyperons (strange baryons), free quarks (forming a hybrid star), or color-superconducting states 
(Alford et al. 2008). It is even possible that the entire star converts into a lower energy self-bound 
state consisting of up, down and strange quarks, known as a strange quark star (Witten 1984; We¬ 
ber 2005), possibly in a solid state (Xu 2003). Other states that have been hypothesized include 
Bose-Einstein condensates of mesons (pions or kaons, the latter containing a strange quark) and 
A baryons (resonant states). The densities at which such phases would appear, and the degree to 
which they might co-exist with other phases, are highly uncertain. 
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In principle, all these phases are governed by the strong force. First principle calculations of 
the interactions for many-body QCD systems are however presently unfeasible due to the fermion 
sign problem. Because fermion wave functions change sign upon exchange of particles, computing 
the expected energy of a system involves the addition of a large number of terms with alternating 
signs. This is especially problematic at high density, and it is at present impossible to obtain direct 
predictions for strongly interacting quark matter. Instead one has to resort to phenomenological 
models. This is particularly true for hadronic many-body systems, since a single hadron is already 
a complex many-body QCD system. Here, the different models are all based on hadronic inter¬ 
actions. For the non-nucleonic degrees of freedom such as hyperons, experimental information 
is scarce and uncertainties even larger than for purely nuclear systems. Lattice approaches and 
effective field theory expansion could help in the future to improve the starting point of ab-initio 
hadronic many-body calculations with these degrees of freedom. Planned laboratory experiments 
on hypemuclei (nuclei containing at least one hyperon) will also give additional constraints. 

Figure 2 compares the parameter space that can be accessed within the laboratory to that 
explored with NS. While the two are complementary, it is obvious that NS provide a truly unique 
exploration space. In addition to the high densities, the gravitational confinement in NS renders 
them sufficiently long-lived such that weak beta equilibrium is achieved and as far as the EOS 
is concerned, NS can be considered as cold. Thus, neither the neutron rich ground state of beta 
equilibrated ultra-dense nuclear matter nor the exotic non-nucleonic states of matter described in 
the previous paragraph can be reached in the laboratory. 

1.2 Unknowns in dense matter superfluidity 

The presence of superfiuid components has significant effects on the star’s long term evolu¬ 
tion and dynamics. Superfluidity suppresses the nuclear reactions that lead to cooling, but also 
leads to new neutrino emission channels. The suggestion of rapid cooling of the young NS in the 
Cassiopeia A supernova remnant (Page et al. 2011; Shternin et al. 2011) has stimulated a detailed 
discussion of the impact on superfluid parameters. If confirmed, such a fast cooling epoch would 
constrain the critical temperature at which the neutrons in the core become superfluid. At the same 
time, the effect requires the onset of proton superconductivity to have taken place much earlier. In 
pulsar timing, superfluidity is thought to be associated with the observed restlessness, ranging from 
the enigmatic glitches seen in young pulsars to the timing noise that appears to be generic in all 
systems. A superfluid rotates by forming an array of quantised vortices (Figure 1). The motion of 
these tiny tornadoes determines the star’s spindown behaviour, e.g. whether it proceeds smoothly 
or in an erratic fashion. 

Key issues involve the friction experienced by the vortices and whether or not they are pinned 
to some other component in the system. Current models suggest that the most important source of 
friction is electron scattering off of magnetic fields associated with the vortices (Alpar et al. 1984; 
Mendell 1991; Andersson et al. 2006). The outcome depends to a large extent on the so-called 
entrainment. This is a non-dissipative effect that determines the dynamical effective mass of the 
particles involved, effectively encoding how mobile the superfluid component is. In the core of the 
star, the entrainment arises due to the strong interaction (Alpar et al. 1984; Borumand et al. 1996), 
while in the inner crust it is due to Bragg scattering off of the nuclear lattice (Charnel 2006). The 
interaction between magnetised vortices and the much more plentiful fluxtubes associated with the 
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Figure 2: The parameter space and the states of matter accessed by NS, as compared to laboratory exper¬ 
iments on Earth. The figure shows temperature T against baryon number density ns against asymmetry 
a = 1 — 2Yq, where Yq is the hadronic charge fraction (in general number of protons/number of baryons). 
a = 0 for matter with equal numbers of neutrons and protons, and a = 1 for pure neutron matter. Green 
region: Isolated nuclei up to a « 0.3, above this one would find a mix of nuclei and light particles. 


proton superconductor (e.g. the star’s bulk magnetic field) may also be dissipative (Link 2003). In 
some regimes it may even cause pinning, which would link the star’s spin evolution to changes in 
the magnetic field. In the star’s crust, vortices are expect to pin to (or between) the lattice nuclei. As 
outlined below, precision timing provides a unique probe of the parameters associated with these 
phenomena, going far beyond the current constraint on the critical temperature from cooling data 
and the ballpark superfluid moment of inertia from glitch observations. 

In the absence of friction and pinning, the neutron superfluid would be more or less decoupled 
from the charged components in the star. This would have a major effect on the moment of inertia 
associated with the spindown, and could in fact lead to a complicated evolution of the braking 
index as the star cools and the superfluid region grows (Ho & Andersson 2012). On the other 
hand, vortex pinning is a key requirement for the “cartoon-level” explanation for glitches, where 
a pinned component builds up an angular momentum reservoir as the system evolves. The stored 
angular momentum is then released in some large scale unpinning event and the pulsar is seen 
to spin up. The relaxation that follows should be determined by vortex friction (Sidery et al. 
2010) and the mechanism that repins the vortices to return the system to the point where it may 
glitch again. This is a very complicated story, and there has been precious little progress in our 
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theoretical understanding of the glitch phenomenon since the first observations at the end of the 
1960s. However, a well-designed SKA campaign could change the situation dramatically. At the 
same time a concerted effort is, obviously, required on the theory side. Finally, superfluidity is 
also expected to have impact on possible free precession. Significant vortex pinning would, in fact, 
prevent slow precession (Jones & Andersson 2001; Link 2003). Instead, one would expect the 
precession period to be similar to the rotation period. Higher quality data may allow us to constrain 
— or even rule out — this possibility. 

2. Methodology 

2.1 Connecting observables to physics using pulsar timing 

The relativistic equations of stellar structure^ relate the EOS of dense matter, which depends 
on the microphysics of composition and the strong force, to macroscopic observables including the 
mass M and radius R of the NS. Figure 3 illustrates the relationship between the EOS and the M-R 
relation for different models: there is a unique map between the two (Eindblom 1992). Constraints 
on or measurements of M and/or R therefore limit EOS models, and hence the interior composition 
of NS. Ideally, we would like to use the SKA to measure M and R for many NS in order to trace 
out the EOS in the M-R diagram. Unfortunately, direct measurements of R in the radio are likely 
impossible, and even indirect measurements through estimates of the moment of inertia are, as 
discussed later in this chapter, incredibly difficult. However measuring masses is something that 
radio observations of pulsars, via high-precision timing measurements, can do exceedingly well, 
and the SKA will be a fantastic tool for this purpose. 

Pulsar timing is the process of unambiguously accounting for every rotation of a NS over 
extended time periods (as long as decades, currently), and it provides astonishingly precise mea¬ 
surements of the pulsar’s spin, astrometric, and (if appropriate) orbital parameters, including how 
each of these changes with time (Eorimer 2008). For “recycled” pulsars, most of which are in 
binary systems, the five Keplerian orbital parameters (orbital period Pj,, projected semi-major axis 
X = ai sin //c (where / is the orbital inclination and c is the speed of light), eccentricity e, time of pe- 
riastron passage To, and the argument of periapsis co) are measured with typically 4—12 significant 
figures, making them essentially perfectly known. 

For some binary pulsars, one or more of the “post-Keplerian” (PK) timing parameters may be 
measurable (Shao et al. 2014). These observables (the precession of the argument of periapsis d), 
the orbital period decay /),, the combination of time dilation and gravitational redshift 7 , and the 
Shapiro delay “range” r and “shape” s) can be expressed in the framework of general relativity (GR) 
as functions of only physical constants (e.g. solar mass in time units Tq = GMq/c^ ~ 5 Jts), the 
extremely well-measured Keplerian parameters, and the masses of the pulsar lUp and the companion 
star mo- Since nip and nic are the only two unknowns, if two PK observables are measured, the two 
equations can be used to solve for nip and nic. Tests of GR come from measurements of additional 
PK parameters. 

*The Tolman-Oppenheimer-Volkoff equations (the equation of mass continuity, force balance, and the EOS P = P(p) 
where p is the total density and P the pressure), modified to take into account rotation. The EOS should also include 
temperature: but for NS, we are so far below Eermi temperature that this can be neglected in computing bulk structure. 
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Figure 3: The link between uncertainty in fundamental nuclear physics, the EOS, and the NS M-R relation. 
The figure shows the pressure density relation (left) and the corresponding mass-radius relation (right) for 
different EOS models. Grey band; range of nucleonic EOS based on chiral effective field theory (Hebeler et 
al. 2013). Black solid: Hybrid EOS from Zdunik & Haensel (2013). Black dashed: Nucleon + hyperon core 
EOS from Bednarek et al. (2012). The other EOS are from Lattimer & Prakash (2001): Red (nucleonic); 
Magenta (quark stars). All models shown are compatible with the existence of 2 Mq NS; discovery of a 
higher mass star would rule out the softer EOS with the lowest maximum masses. 


Unfortunately, only a small fraction of binary pulsars have timing precisions, orbital parame¬ 
ters, or companions which allow the measurement of PK parameters. To measure 6), for instance, 
which by itself provides the total system mass mp + nic and is exactly the same effect as Mercury’s 
precession of perihelion, we need pulsars in eccentric enough orbits so that we can trace ft) pre¬ 
cisely. Most millisecond pulsars (MSPs) in the Galactic disk have extremely circular orbits and so 
even though their orbital parameters are measured to extraordinary precision, ft) cannot be. How¬ 
ever, globular clusters contain many eccentric MSPs (Ransom 2008; Hessels et al. 2014) and there 
is evidence for a population of them in the Galactic disk as well (Champion et al. 2008; Barr et al. 
2013), providing opportunities for new discoveries of eccentric systems by the SKA. 

Detection of orbital period decay Pf,, due to the emission of gravitational radiation, demands 
compact systems such as the original Hulse-Taylor binary (Weisberg et al. 2010). These double 
NS (DNS) systems (Kramer et al. 2006; Fonseca et al. 2014) tend to be in highly eccentric orbits, 
which allows measurement of ft), and after enough precession has occured, y (the requirement for 
measurable precession makes y effectively impossible to measure in systems with near-circular 
orbits, even if they are compact and relativistic). Several compact NS—white dwarf (WD) systems 
have been discovered recently which allow measurement of (Bhat et al. 2008; Freire et al. 2012; 
Antoniadis et al. 2013), implying that the SKA should find many more such systems. 

The final PK parameters of interest are those making up the Shapiro delay, which is a rel¬ 
ativistic propagation delay, typically of the order of /fs, caused as the pulses travel through the 
gravitational potential of the companion star. The amplitude of this effect depends sensitively on 
the inclination angle and requires nearly edge on viewing and high timing precision to become 
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measurable. With huge improvements in timing precision over the last decade, as well as the 
discovery of many new MSPs, the number of Shapiro delay measurements has increased substan¬ 
tially, and that number should increase tremendously with the SKA. The first high-precision 2Mq 
NS was measured with this technique (Demorest et al. 2010), and that single number has provided 
exceptionally strong constraints on the EOS (Hebeler et al. 2013). 

For non-binary slowly-spinning pulsars it is impossible to measure orbital or PK parameters, 
yet rapid-cadence (daily or every few days) timing observations can detect the timing “glitches” 
which young pulsars often exhibit (Section 3.4). These virtually instantaneous increases in the 
pulsar’s spin rate and (sometimes) spin-down rate, and the subsequent relaxation process, probe 
the superfluid interiors of NS. 

2.2 Importance of large surveys to sample populations 

Surveys continue to be a major driver for pulsar science and its applications. By increasing the 
known population, we increase the number of outlier systems whose exceptional properties make 
them powerful physics laboratories, e.g.(Kramer et al. 2006; Ransom et al. 2014). Consider the fact 
that there are 2300 known radio pulsars (Manchester et al. 2005) out of an estimated total Galactic 
population of a few tens of thousand sources beaming towards the Earth (Keane et al. 2014). Of 
these known sources, fewer than 10 % are in binary (or trinary) systems, a pre-requisite for preci¬ 
sion mass measurement through either Shapiro delay or a combination of PK orbital parameters 
and constraints from spectroscopic modeling of the optical companion star. Unfortunately, but un¬ 
surprisingly, only a handful of systems have extreme properties that push the envelope in terms of 
constraining the EOS. There are only 10 radio pulsars with rotational periods less than 2ms, and 
only 2 NS with precision mass measurements of Mns ^ 2 M 0 (Demorest et al. 2010; Antoniadis 
et al. 2013). Pushing to faster spin rates, which can constrain the maximum NS radius, and higher 
masses will require both a significant increase in the known population (a doubling or tripling) as 
well as careful follow-up and precision timing of the most promising discoveries. 

In addition to strong constraints on the EOS from maximum masses, a much larger sample 
will also provide precision mass measurements across the known distribution Mns ~ (1.2 — 2)Mo 
(lower mass limit from Janssen et al. (2008); Ferdman et al. (2014)). Mapping the NS mass dis¬ 
tribution (Tauris et al. 2014) is important for understanding the a^^trophysics of NS formation in 
core-collapse supernovae (EOS dependent) and subsequent ‘recycling’ in a binary system (Alpar 
et al. 1982). The lowest mass NS are also potentially extremely interesting. If these systems are 
formed by electron capture supernovae (Kitaura et al. 2006; Ferdman et al. 2014), this could place 
a tight constraint on the EOS in terms of the relation between gravitational and baryonic mass 
(Podsiadlowski et al. 2005). This idea is, however, subject to uncertainties in both the astrophys- 
ical setup (such as the progenitors), and on the nuclear physics side with respect to the nuclear 
reactions during the SN and the EOS used for the extraction of the masses. 

SKAl-MID and SKAl-EOW will conduct both wide-field and targeted pulsar surveys. It is 
expected that an all-sky survey will discover ~ 10000 pulsars (Keane et al. 2014), of which ~ 1800 
will be MSPs and ~ 100 DNS. Though impossible to predict precisely, ~ 10% of these should be 
well-suited for EOS constraints. SKAl will also do deep targeted searches of unidentified Fermi 
7 -ray sources (Keane et al. 2014) and the Galactic globular cluster systems (Hessels et al. 2014). 
These searches will discover primarily MSPs, some of which will also be EOS laboratories. The 
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Figure 4: Current NS mass measurements, 68% confidence intervals (http;//jantoniadis.wordpress.com/ 
research/ns-masses/ for original references). Blue -LMXBs; Red - Eclipsing MSPs; Green - DNS; Magenta 
- NS-WD binaries; Black - MSP-Main Sequence binaries; Cyan - Triple Pulsar. Though some sources may 
have higher masses than PSRs J1614-2230 and J0348H-0432, there are large systematic uncertainties on these 
measurements or other reasons why they are not as direct or reliable a measurement. 

current spin record holder resides in the globular cluster Terzan 5 (Hessels et al. 2006). A signifi¬ 
cant fraction of the MSPs found in targeted searches will be in eclipsing binaries (‘black widows’ 
and ‘redbacks’ (Roberts 2011)), which, although more difficult to detect and time precisely, may 
harbour the most massive NS created by nature (van Kerkwijk et al. 2011; Romani et al. 2012) be¬ 
cause they have accreted significant mass from their companion and/or were formed more massive. 

3. State of the art (observations and theory) 

3.1 State of the art: masses 

A measurement of a NS mass, even without a simultaneous measurement of R, constrains the 
EOS if it exceeds the maximum mass predicted by a given model (Figure 3). Figure 4 shows current 
NS mass measurements. The two most constraining systems are PSRs J1614—2230 (Demorest 
et al. 2010) and J0348-I-0432 (Antoniadis et al. 2013), with masses of 1.97 ±0.04 and 2.01 ± 
0.03 Mq respectively. PSR J1614—2230 is a 3.1 ms MSP orbiting an intermediate-mass WD every 
8.7 days. The system has a very high inclination (89.17°) resulting in a strong Shapiro delay signal. 
PSR J0348-I-0432 is a 39 ms pulsar in a relativistic, 2.46 h orbit with a low-mass WD. For this binary 
the PK parameters cannot be measured and the masses have been measured through phase-resolved 
optical spectroscopy of the bright WD companion. 

Current knowledge of hadronic physics suggests it is nearly impossible for matter to be purely 
nucleonic at the expected massive NS core densities of (7-10) Psat- No matter what the precise 
composition, all alternative possibilities add additional degrees of freedom. The immediate effect is 
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a softening of the EOS, lowering the maximum mass, unless the interaetion of the new eonstituents 
provides the neeessary repulsion to eounter-balanee the strong gravitational foree. While 2Mq NS 
do not direetly exelude any eomponents in the NS interior they do, however, put very stringent 
eonstraints on the respeetive interaetions. In partieular, they have triggered intensive diseussions 
on the "hyperon puzzle": ab-initio many-body ealeulations prediet the appearanee of hyperons at 
2-3 Psat but the maximum masses of sueh models do not exeeed ~ 1.7 M©. Proposed solutions 
inelude a very early transition to quark matter, additional repulsion in the ab-initio ealeulations via 
two- or more-body hyperonie forees and the appearanee of other new partieles sueh as A-baryons. 

Evidenee for even more massive NS has been reported for eelipsing binary pulsars with low- 
mass eompanions (van Kerkwijk et al. 2011; Romani et al. 2012). Based on optieal observations 
of the eompanions, these systems seem to host NS with M ~ 2.4M0, well above the predietions 
of many EOS models. However, there are fairly large systematie uneertainties assoeiated with 
lighteurve modeling (whieh mostly impaet the determination of inelination), rendering these mea¬ 
surements so far unreliable for EOS eonstraints. A signifieant improvement of these measurements 
will most likely require synergy of sensitive radio and optieal observations with next-generation 
observatories sueh as SKA and E-EET (Antoniadis et al. 2014). 

3.2 State of the art: moments of inertia 

Moment of inertia / is a funetion of M and R, so yields R if M is known independently. 
Current eonstraints on I are very approximate, and are based mostly on y-ray flux measurements. 
PSR J1614—2230 is an espeeially interesting ease, sinee its mass is well known. This pulsar was 
among the first MSPs deteeted in y-rays (Abdo et al. 2009). Assume that the spin down luminosity 
E = An^lP/P^ (Eorimer & Kramer 2005). Sinee the measured y-ray luminosity must be < 100% 
of E, this yields a lower limit on 7. If the distanee to PSR J1614—2230 is given by the free-eleetron 
model of the Galaxy (Cordes & Eazio 2002), this yields 7 > 10"*^gem^ (not very eonstraining). 
Distanee estimates, whieh fix the y-ray luminosity, are however uneertain. The distanee is also 
neeessary to evaluate kinematie eontributions (whieh may be large) to P. Detailed knowledge and 
modelling of the emission geometry is also neeessary to estimate the real luminosity of the pulsar: 
the estimate above assumes isotropie emission, whieh is not appropriate sinee in this ease there 
would be no pulsations. One advantage of this method, however, is that sueh lower limits on 7 
will eventually be available for a variety of other pulsars (Guillemot et al. 2012) with a variety of 
masses, something that might improve its stringeney. 

3.3 State of the art: spin rates 

During reeyeling, a bare non-magnetie NS eould easily attain very rapid spin after aeereting 
only ~ O.IM© (Cook et al. 1994; Burderi et al. 1999)^ . Irrespeetive of the proposed EOS, the 
spin frequeney / must be lower than the Keplerian frequeney fx above whieh a NS, due to een- 
trifugal forees, beeomes unstable to mass shedding at its equator. To a very good approximation 
ffc = C(M/M0)^/^(7?/lOkm)^^/^ (Haensel et al. 2009). where M is the (gravitational) mass of the 
Keplerian eonfiguration and R denotes the (eireumferential) radius of the non-rotating eonfigura- 
tion of the same gravitational mass. C is 1.08 kHz for NS and 1.12 kHz for strange stars (Haensel 

^Although there are limiting mechanisms that might halt spin up, see Patruno & Watts (2012) for a review. 
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Figure 5: Left: Constraints arising from different spin rate measurements. NS of a given spin rate must 
lie to the left of the relevant limiting line in the M-R plane (shown in blue for various spins). EOS models 
as in Figure 3. The current record holder, which spins at 716 Hz (Hessels et al. 2006) is not constraining, 
but given a high enough spin individual EOS can be ruled out. Between 1 kHz and 1.25 kHz, for example, 
individual EOS in the grey band of nucleonic EOS would be excluded. Right: The spin distribution of radio 
MSPs with a spin frequency > 100 Hz. The spin distribution of the smaller sample of accreting NS is similar 
(Patruno & Watts 2012; Watts 2012). 


et al. 2009), obtained from fits to general relativistic calculations of rotating stars. Since we require 
/ < Ik, this results in a constraint on the radius: 


R<\0 


■ M ■ 

1/3 

[ / 1 



1 kHz 


1 -2/3 


km 


(3.1) 


The more compact the star is, the higher the supported rotation rate can be, and the tighter the 
constraint on the EOS. Figure 5 shows the spin distribution of the known radio pulsars. The binary 
MSP J1748—2446ad (Hessels et al. 2006) is the NS with the shortest known rotational period, 
Pmin = 1-396 ms, corresponding to a spin frequency of 716 Hz (the most rapidly rotating accretion- 
powered NS rotates at 619 Hz (Hartman et al. 2003)). The resulting constraint is shown in Figure 
5: the radius of a non-rotating 1.4 Mq star would have to be smaller than 15 km. This is not 
particularly constraining for current models: rotation rates faster than a millisecond are required to 
rule out EOS. 


3.4 State of the art: glitches 

The majority of the glitches known today have been detected in the data of long monitoring 
programs (at Jodrell Bank, Par'kes and Urumqi; (Yuan et al. 2010b; Espinoza et al. 2011; Yu et al. 
2013)). This is because glitches are inferred from observations before and after the glitch epoch. 
The most minimal characterisation of a glitch is the measurement of the step in spin frequency (Av). 
However, many glitches also induce a change in spin-down rate (Av) and a process of relaxation 
towards the pre-glitch rotational values, with timescales that can go up to a few hundred days. 
Constraints and information on the glitch mechanism and the rotational dynamics of the neutron 
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Figure 6: Number of glitches per year for most pulsars known to have glitched and observed for more than 
3 yr (by 2010, from Espinoza et al. (2011)). The straight line is a linear ht to the maximum value of Ng in 
each half-decade of characteristic age Tf. 


superfluid can be obtained if we can detect and characterise many (if not all) glitches in several 
sources, precisely. Our ability to do this depends directly on the cadence of the observations, 
as well as on the accuracy of the TOAs. At the moment we know of ~ 440 glitches in ~150 
pulsars. Glitch activity is concentrated in young objects (Tc < 100 kyr) and there is an almost linear 
correlation with spin-down rate, i.e. the faster a pulsar spins down, the higher its glitch activity 
(Lyne et al. 2000; Espinoza et al. 2011). The occurrence of glitches is in general not periodic and 
most pulsars exhibit less than one glitch every 2 years. This shows the importance of monitoring a 
large number of pulsars over a long time (Figure 6) using sub-arraying to observe multiple sources 
at higher cadence when possible. 

The size of the largest glitches has already constrained the amount and location of the partic¬ 
ipating superfluid (Andersson et al. 2012; Charnel 2013; Piekarewicz et al. 2014). The detection 
of small glitches and the lower end of the Av distribution are however highly uncertain. Below 
certain sizes detection is compromised and there is confusion with phenomena like timing noise 
and disperson measure (DM) variations. A study of the glitch activity of the Crab pulsar has shown 
that the number of glitches decreases significantly below Av = 0.05 /f Hz, implying a rather large 
minimum glitch size (Espinoza et al. 2014) which can in principle constrain the trigger mecha¬ 
nism. However, the minimum size found for the Crab glitches is larger than many glitches detected 
in other pulsars. In-depth studies of other pulsars will be important to assess the universality of the 
existence of a lower cut-off and to reach a better comprehension of the glitch trigger mechanism. 
To accomplish this it is necessary to have observations sensitive to very small glitches: the Crab 
pulsar glitch study was enabled by the high cadence (one observation per day) and high precision 
of the data. 

Glitch recoveries provide information on the dynamics of the internal superfluid and its in¬ 
teraction with the other NS constituents. The Vela pulsar is the prototype source showing glitch 
recoveries, which are also observed in many other young pulsars. During the relaxation process the 
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Figure 7: Observable Shapiro delay amplitude - i.e., the portion of the signal that is not covariant with other 
orbital parameters - as a function of inclination i and binary companion mass. Assuming orientations have 
uniform cosi, at timing precision of ~50 ns, >80% of NS-WD binary systems with nic >0.1 M© should 
show detectable Shapiro delay. 

spin-down rate slowly decays from the rather large (very negative) value induced by the glitch. It 
takes years for it to again reach pre-glitch levels and in many cases the process will be interrupted 
by a new glitch. PSR B2334-I-61 exhibited one of the largest glitches known to date, followed by 
a standard recovery process. Nonetheless, significant oscillations were reported in post-glitch tim¬ 
ing residuals (Yuan et al. 2010a). To date no post-glitch oscillations have been detected in other 
pulsars. This could be related to the size of the glitch but could also just be an observational bias 
against detecting similar effects in lower sensitivity data taken at lower observing cadence.. 

Other pulsars have shown peculiar glitch recoveries, like the RRAT PSR J1819—1458 and 
the radio pulsar PSR J1119—6127 (Lyne et al. 2009; Weltevrede et al. 2011). In these two cases 
the spin-down rate over relaxed and reached values significantly smaller than the pre-glitch lev¬ 
els. Pulse shape variations were detected in a couple of observations following the glitch in PSR 
J1119—6127 but the rather low cadence available (month^^) prevented a more detailed study. The 
unique recoveries were attributed to the large dipole fields of fhe fwo pulsars, fhough no physical 
mechanism was proposed. The analysis of glifch recoveries can fell us abouf differences befween 
fhe regions responsible for fhe rofafional relaxafion, for insfance whefher fhere is pinning in fhe 
oufer core or nof. These sfudies would benefif from higher cadence and more precise observations. 

4. Expectations from SKA 

4.1 Improvements in mass measurements 

The SKA is expecfed fo improve on fhe number of measured NS masses in fwo ways. Firsf, 
fhe dramatic increase in fhe number of known radio pulsars will yield more exotic sysfems in which 
masses can be defermined (e.g. via measurable Shapiro delay). We mighf expecf fhese fo comprise 
roughly fhe same fracfion of fhe fofal populafion as is observed now, or ~ 1-2%. Even so, fhis effecf 
alone would increase fhe number of fofal mass measuremenfs by a facfor ~ 10. The second way 
fhaf fhe SKA will improve mass measuremenfs is fhrough fhe improved timing precision fhaf will 
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be achievable with the massive increase in telescope sensitivity. This will improve the uncertainty 
on already-measured masses. It also has potential to increase the fraction of systems in which 
masses can be measured, boosting the total numbers even further. Detection of Shapiro delay in 
significant fractions of the NS-WD population may become possible (Figure 7). 

In contrast to many of the other relativistic timing parameters, Shapiro delay is not heavily de¬ 
pendent on special environmental or evolutionary history. The effect applies to all binary systems, 
is insensitive to contamination by non-relativistic effects, and its detectability mainly depends on 
random geometrical orientation - at high inclinations, the radio pulses we receive must travel close 
to the companion star, resulting in a much larger time delay^. Furthermore at low inclinations, the 
effect becomes highly covariant with the Keplerian orbital parameters, which reduces its observ¬ 
able amplitude even further (Freire & Wex 2010). Nevertheless, if high enough timing precision 
can be achieved, the effect can be detected even in fairly low-inclination systems. This is illus¬ 
trated in Figure 7, which shows the observable (non-covariant) Shapiro delay amplitude versus 
cos/. At ~50 ns timing precision, Shapiro delay is potentially detectable at inclinations as low as 
~ 40° - this includes ~80% of all binary systems assuming random orientations (uniformly dis¬ 
tributed in cos/). The SKA is expected to achieve sufficient timing precision, and discover enough 
new binary systems, that we might expect at least a factor of five more Shapiro delay based mass 
measurements. 

4.2 Improvements in moment of inertia measurements 

The SKA opens up the possibility of a precise measurement of moment of inertia I for one 
pulsar for which the mass is already precisely known, the Double Pulsar PSR J0737—3039A. There 
are three routes by which pulsar timing measurements can lead to a measurement of I (Kramer & 
Wex 2009). A highly relativistic system such as the Double Pulsar (Burgay et al. 2003; Lyne et 
al. 2004) is required in order for the effects of I on the pulsar timing to be discernable. The first 
measurable effect is expected to be a contribution comparable in size to the 2PN-level correction 
to the advance of periastron, 6). This would produce an offset between the measured value of d) 
and that expected by taking the other measured PK timing parameters and combining them for 
a (IPN) prediction. In practice, this will require very precise measurement of at least two other 
PK pai'ameters which produce largely orthogonal curves in the mass-mass plane. For the Double 
Pulsar, the parameters best suited to this prediction will be the sine of the orbital inclination angle, 
s, and the rate of decay of the orbital period, Pb. Simulations assuming 5 ps timing precision predict 
a 10% measurement of I in 20 years (Kr'amer & Wex 2009). Timing with SKAl should do even 
better than this. An important point is that the corrections to Pb for the relative acceleration of the 
Solar System and the pulsar' system must be known to high accuracy in order for this method to 
work. These corrections require precise measurements of the pulsar’s proper motion and distance, 
things that should be easily achieved for the Double Pulsar with the timing and VLBI capabalities 
of SKAl. A further issue is the need for a reliable model of the Galactic potential; as the Double 
Pulsar is close to the Plane, this should not pose an insurmountable problem in this particular case, 
and models will also be improved on the necessary timescale by the results of the GAIA satellite. 

^Note however that detectability of a Shapiro delay would be adversely affected if the inclination were so high that 
the system was in fact eclipsing. 
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The second method involves determining non-linear time evolution of the longitude of perias- 
tron, if the spin of the (recycled) pulsar is inclined with respect to the system’s angular momentum. 
This should allow the measurement of the projection of the pulsar’s spin vector on the orbital plane, 
leading to I via a precession-determined estimate of the spin-orbit misalignment angle 5. The third 
technique relies on measuring two time derivatives of the projected semi-major axis of the pulsar’s 
orbit, V, once again leading to an estimate of 7 sin 5, although this effect would need to be sepa¬ 
rated from aberration-induced changes in projection of the orbit (Damour & Schafer 1988; Wex & 
Kopeikin 1999; Damour & Taylor 1991). For the Double Pulsar these last two techniques are not 
promising, as the spin-orbit misalignment angle for the recycled pulsar is known to be small (Ferd- 
man et al. 2013) and furthermore the orbit happens to be nearly edge-on to the line of sight, making 
the changes in projection tiny. The small misalignment angle is likely the result of a fairly symmet¬ 
ric, low-mass-loss explosion in the second supernova in the system (Podsiadlowski et al. 2004), but 
there are other DNS systems that are known to be significantly misaligned (Kramer 1998; Weis- 
berg & Taylor 2002; Fonseca et al. 2014). The SKA (1,2) is expected to discover (100,180) DNS 
systems, and it is reasonable to hope that favourable geometries, and indeed even more relativistic 
orbits will be found in some of them. 

4.3 Improvements to spin measurements 

In the more than three decades since the discovery of the first MSP B1937-1-21 {P = 1.56 ms) 
(Backer et al. 1982), only one object, PSR J1748—2446ad {P = 1.39 ms) (Hessels et al. 2006) has 
been found that spins faster. Until a few years ago, the sometimes prohibitive computational costs 
of a large survey sensitive to very rapidly rotating objects lent support to the hypothesis that there 
might have been an observational bias limiting the number of known MSPs with spin period close 
to that of the record holders or below. However, the improvement of the computational capabilities 
implemented in modern pulsar search experiments made them nominally able to pick up ultra-fast 
radio pulsations and several new binary MSPs with spin periods below 2 ms have been discovered 
in recent years. An additional contribution to the increment of this population of rapidly spinning 
NS came from the targeted radio searches towards y-ray unidentified sources selecfed from fhe 
Fermi cafalogues. 

Inferesfingly, in 9 cases (over a fofal of 11 binary MSPs having P <2 ms), including fhe 
spin record-holder PSR J1748—2446ad, fhe mosf rapidly rofafing pulsars appear fo be hosted in 
eclipsing sysfems: binaries in which fhe mailer from fhe companion obscures fhe radio pulsalions 
for a large fraclion of or perhaps even fhe entire orbil. This suggesls lhal allhough Ihere may be 
no observational bias againsl discovering fhe mosf rapidly relating systems, Ihey be nonelheless 
be hidden by obscuration from Ihe companion slar. However, Ihe eclipses in Ihese systems appear 
lo be slrongly frequency dependenl, wilh Ihe eclipses being signilicanlly longer al frequencies al a 
GHz and below. Il is also interesting lo note lhal many pulsars in eclipsing systems were discovered 
in a 2 GHz survey of globular clusters wilh Ihe GBT. 

Due lo Ihe very large numbers of beams needed lo cover Ihe sky and accounting for Ihe long 
dwell time needed lo reach Ihe required sensitivity (a pulsar like J1748—2446ad has a flux den¬ 
sity of only ~ 80 /tJy al 2 GHz) a deep all-sky survey al 2 GHz is nol feasible wilh Ihe currenl 
generation of radio telescopes. The beam-forming capabilities (allowing simullaneous search of 
tens or hundreds of sky positions), and Ihe higher inslanlaneous sensitivity of SKAl-MID (reach- 
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ing ~ 10 /^Jy sensitivity in a relatively short integration time) will give the first opportunity for a 
deep all-sky survey at 2 GHz"*. Provided the sub-millisecond pulsars exist and are mostly hidden in 
eclipsing binaries, this experiment would be the best designed ever in the radio band with the aim 
of uncovering at least the brightest objects in this category. A similar later experiment performed 
with SKA2 will be finally sensitive even to sub-millisecond pulsars in a portion of the Galaxy well 
beyond the volume probed by the Fermi y-ray satellite. 

Even if SKAI-MID and SKA2 do not find sub-millisecond pulsars, SKAI-MID will certainly 
discover at least few MSPs spinning faster than PSR J1748—2446ad, while SKA2 will deliver a full 
census of the population of the of the nearby, thus unaffected by scattering and dispersion delays, 
ultra rapidly spinning MSPs. Likely, the accurate timing observations possible with SKA will also 
lead to a measurement of the pulsar mass (via relativistic effects) for at least a sub-group of the 
objects. It will then be possible to exploit the large sample of rapidly spinning sources in order 
to derive statistically sound results. Population analysis (Possenti et al. 1999; Chakrabarty 2008; 
Ho et al. 2011; Papitto et al. 2014) may constrain not only the EOS, but will also elucidate other 
fundamental physical processes related to the final stage of the recycling that may depend on the 
EOS e.g.: (i) the role of the emission of gravitational wave radiation (r-modes or accretion induced 
quadrupole moments (Bildsten 1998; Bondarescu & Wasserman 2013)) in halting the pulsar spin 
up; (ii) the physical process provoking the decay of the surface magnetic field in the MSPs; (Hi) 
the relative importance of the radio ejection (Burderi et al. 2001) and of the Roche-lobe decoupling 
phase (Tauris 2012) in shaping the spin distribution of the MSPs (see also Kiziltan & Thorsett 
( 2010 )). 

The impact of a NS rotating at over ~ 1.1 kHz would be substantial, and 1.5 kHz would be 
revolutionary: no current models can sustain such a fast spin. A hypothetically observed frequency 
of 1.4 kHz, for example, would constrain the radius of a 1.4 Mq NS to less than 9.5 km (Eigure 5). 
In general higher frequencies favour additional non-nucleonic degrees of freedom, since they give 
more compact stars. Simultaneous measurement of high mass and high spin would be particularly 
constraining (Haensel et al. 2008). While massive NS require repulsion among their constituents 
at high density, fast rotation requires the interaction not to be too repulsive to counter-balance the 
centrifugal force. 

4.4 Improvements in glitch sampling 

The detection of glitches and the characterisation of their recoveries rely strongly on the con¬ 
stant monitoring of a sample of pulsars. The observing capabilities which are currently offered 
by the radio observatories that monitor young pulsars are sufficient for the detection of medium 
to large glitches. The detection of small glitches and the study of the lower end of the glitch size 
distribution are, however hard to attain with the current observing programs. SKA glitch studies 
should focus on monitoring a selected group of pulsars and should be able to provide the com¬ 
plete glitch sample for the observed sources during the observing campaign. Besides uncovering 
the lower end of the size distribution, this will allow the study of waiting time (inter-glitch time 
intervals) distributions, which can also shed light on the glitch mechanism. 

"^If observing time does not permit a full survey, one would instead focus on pointed observations of GCs, Fermi 
sources, and quiescent LMXBs to find the fastest spinning pulsars. 
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There are two main factors limiting the detection of glitches (Espinoza et al. 2014): 1) How 
often pulsars are observed, and 2) o^, the larger of the dispersion of the timing residuals in a 
given time-interval (10-15 observations) and the average uncertainty of the TOAs. Compared to the 
current observing setups the SKA will offer superior sensitivity, hence we assume that will never 
be dominated by TOA uncertainties. Unmodelled rotational irregularities and observing cadence 
will likely be the main limiting factors in glitch detection. Among the known effects that can 
increase the dispersion of the timing residuals are DM variations and timing noise. DM variations 
will be corrected via regular multi-frequency observations (highlighting the need for SKAl-LOW), 
necessary to achieve high precision timing, and therefore should not affect glitch detection. Timing 
noise, on the other hand, still needs to be well understood in order to correct the data for its effects. 
The detection of small glitches (for a given pulsar) will therefore be dependent on the intrinsic 
levels of timing noise and its characteristic timescales in combination with the observing cadence. 

The nominal cadence of SKA observations for new pulsars is 1 observation per 2 weeks, for 
at least a year. In addition, timing campaigns for pulsars should last at least 10 years. Although 
this might seems appropriate, with one observation every two weeks the SKA might not be able to 
improve substantially the statistics of small glitches already accomplished by some of the current 
programs, nor to sample satisfactorily post-glitch recoveries. Glitch studies need a higher cadence 
and a longer campaign over a fixed set of pulsars. More frequent SKA observations are planned 
for a sub-sample of interesting MSPs, but these are expected to glitch at an average rate lower than 
1 glitch per millennium (Espinoza et al. 2011). More complete monitoring of a larger sample of 
pulsars may be achievable using SKAl-EOW, which, with the wider field of view and sufficienlly 
large numbers of fied-array beams, should facililale daily observafions of large numbers of sources. 

Glifch defeclabilify limifs also depend on fhe size of fhe spin-down rale sfeps, improving for 
smaller steps. Known glilches have |Av| values in fhe range 10^^ lo 10"^ (x 10 ^^Hzs ^). Thelwo- 
weeks cadence sfrongly limifs fhe defeclion of glilches wifh |Av| > 1 x 10^^^ Hz s^^ Eor insfance, 
some of fhe Crab pulsar glilches could nol be unambiguously delecled, as demonsfrafed in Eigure 8. 
The bias in glifch defeclion is clear from fhe leff hand side plof, being dominafed by al low | Av| 
values. Eor some pulsars is driven by liming noise, buf for ofhers is driven by TOA precision. 
The righf hand side plof shows fhe increased paramefer space lhaf SKA observafions could achieve 
wifh daily observafions and a liming precision of al leas! 5 x 10^^ rolalions (equivalenf lo 5 x 
10^ ns, for a period of 0.1 s). The SKA could easily achieve fhis precision for many pulsars, buf fhe 
level of liming noise presenl in a given fime-inlerval (10-15 TOAs) could be higher. 

The number of pulsars monifored for glilches need no! be large (20-30 objecfs) and should 
include fhose lhaf hav e already been monifored for years (lo Irack long-term glifch hisfory), which 
are brighl for SKA slandards. Precision can Ihus be Iraded for higher cadence. The Aperlure 
Array syslems should allow an optimum combinalion of sensilivily, field-of-view and number of 
beams lo yield exceplional cadence on a very large number of sources. Wifh one short (~ 2 mins^) 
observalion per day giving a TOA precision of al leasl 5 x 10^^ rolalions, over years and for a 
sample of pulsars, we can explore fully Ihe glifch size dislribulion and Irack Ihe recoveries of all 
detected glilches. SKA will Ihus yield Ihe ideal sample lo probe Ihe glilch mechanism. 

^Assuming A^ff/Eys = lOOOm^ (SKA 1) and for a duty cycle of 10%, average flux 0.5 mjy and a bandwidth 

of 500 MHz, 2 mins, of observation gives 2 x 10^^ s precision (for P=0.1 s) 
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Figure 8: 245 glitches in the |Av|-Av plane (Espinoza et al. 2011). Open diamonds: Crab pulsar glitches. 
Left: Detectability limits for observing cadences of 14 days (solid thick line) and 30 days (solid line). 
The segmented lines are for equal to 0.01 (rightmost) and 0.001 (leftmost) rotations, which are values 
commonly reached today. Right: Same as Left but for a cadence of 1 day and TOA precision of 5 x 10^^ 
rotations. 

5. Summary 

The nature of the strong interaetion is fundamental to understanding the behaviour of sys¬ 
tems from nuelei to neutron stars, and remains a serious ehallenge to physieists even though it is 
mathematieally well-defined. SKAl, and ultimately the full SKA, will progressively give us the 
opportunity to explore the nature of matter, the strong foree and superfluidity under the extreme 
eonditions that prevail inside NS, via preeise measurement of pulsar masses and moments of iner¬ 
tia, the diseovery of rapidly-spinning and new types of pulsars, and high eadenee observations to 
deteet pulsar glitehes. 
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